Analytical treatment of 2D steady flames anchored in high-velocity streams 
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The problem of burning of high- velocity gas streams in channels is revisited. Previous treatments 
of this issue are found to be incomplete. It is shown that despite relative smallness of the transversal 
gas velocity, it plays crucial role in determining flame structure. In particular, it is necessary in 
formulating boundary conditions near the flame anchor, and for the proper account of the flame 
propagation law. Using the on-shell description of steady anchored flames, a consistent solution 
of the problem is given. Equations for the flame front position and gas- velocity at the front are 
obtained. It is demonstrated that they reduce to a second-order differential equation for the front 
position. Numerical solutions of the derived equations are found. 

' PACS numbers: 

<D '. 

V/-) \ I. INTRODUCTION 

(N ; 

The problem of flame propagation in a flow of gaseous mixture with a fixed ignition point is important from both 
experimental and theoretical points of view. In fact, relative simplicity of practical realization and importance in 
applications made anchored flames one of the most popular topics in combustion science. Despite these circumstances 
theoretical description of the subject is far from being complete. It would not even be exaggeration to say that 
the very mechanism of formation of steady flame configurations is not fully understood. One of the most difficult 
problems here is the influence of anchoring system on the flame, in particular, the question of its locality. Analytical 
investigation of this arduous question is so complicated that it is usually not raised at all. Another closely related 
problem is the identification of mechanisms driving the development of flame disturbances. It is known for a long time 
that anchored flames may develop bulbous structures their stability properties are affected by gravity diversely 
and deeper than in the case of freely propagating flames 0, H[ , which means that the stabilizing mechanisms in the 
^ j— I * two instances are quite different. 

These issues are especially nontrivial in the case of two-dimensional flames. Namely, reduction of dimensionality 
l— ~~ '■ changes the long-range behavior of the Green functions, leading to enhanced interaction between distant parts of the 
I | flame. Specifically, pressure distribution is determined by the Green function of the Laplace operator, whose integral 
. kernel grows logarithmically with distance in two dimensions, and therefore so does the response to a point source. 
^\ ' This indicates that the regions with large velocity gradients, in particular, vicinity of the anchor, may have strong 
, nonlocal impact on the global flame structure. 

Of particular interest are the high-velocity streams. Experiments indicate that when the incoming gas velocity 
significantly exceeds the normal flame burning speed, the flame front assumes highly elongated shape which often 
can be well approximated by straight lines (V-flames). Even if the front shape is not piecewise linear, simplifications 
admitted by the high- velocity limit make it accessible for theoretical investigation. An important example is the flame 
anchored in a high-velocity uniform stream in a channel. The first detailed theoretical account of this case was given 
by [H who derived an integral equation for the front position, and obtained its numerical solutions [the main results 
. of this work are reproduced in the book The problem was considered independently by Q, whose results were 
J^j ' subsequently critically reviewed and clarified by In the latter work, in particular, the main assumptions employed 
rS I in the analysis were identified and used to derive an integral equation similar to that of 
■ All three works deal with the steady case and assume the following: 

1) For sufficiently high stream velocity, curvature of the stream lines can be neglected. The gas velocity is thus 
parallel to the channel walls everywhere. This implies that the gas pressure is constant in every cross section of the 
fresh and burnt gas regions. Furthermore, neglecting the relatively small constant pressure jump across the front 
makes pressure constant in every cross section of the channel. 

2) The fresh gas velocity is also constant in every cross section, both up- and downstream of the ignition point. 
Neglecting velocity jump at the front makes the flow field continuous everywhere in the channel. 

Using the Bernoulli integral under these assumptions, it is straightforward to show that the pressure is monotonic 
along the channel, so that it can be taken as an independent coordinate, [seeQ. Alternatively, as an independent 
variable can be taken the fresh gas velocity which is also monotonic along the channel, [see 0] • In either way, the use 
of mass conservation gives an integral equation for one of the flow variables. 

It should be mentioned that no attempt was made in the cited papers to justify the approximations 1), 2) any more 
rigorously than outlined above. Although these approximations look naturally, exclusion of the transversal velocity 
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component from the list of dynamical variables represents a quite nontrivial step. Indeed, replacing the system of two 
Euler equations by the single Bernoulli integral means that the two velocity components are completely decoupled 
from each other. Considered on its own, this reduction of the system is legitimate under the assumption of high 
stream velocity. But the flow equations themselves do not constitute the complete system of equations governing 
flame propagation. They must be supplemented by the evolution equation and the jump conditions at the flame 
front, as well as boundary conditions at the channel walls. But the two velocity components are strongly coupled 
by the evolution equation: relatively small variations of the transversal component give rise to large variations of 
the longitudinal component [Cf. Eq. (|33|) below]. At the same time, exclusion of the transversal component from 
the consideration based on the assumptions 1), 2) makes this essential equation obsolete. Thus, in order that the 
evolution equation take its proper place in the analysis of flame propagation, it is necessary to bring the transversal 
velocity component back into consideration. This requires in turn restoration of the remaining flow equation and the 
corresponding jump condition. 

An important step in validating the model described above was made by [8| who showed that the equations derived 
by Zel'dovich and Tsien are asymptotically exact. More precisely, he formulated the way the high- velocity limit should 
be taken in the complete system of governing equations, and demonstrated that the resulting system reduces to an 
integral equation which is equivalent to the equations derived by these authors. One of the purposes of the present 
paper will be to show that Cherny's consideration is not complete. Namely, it omits one of the boundary conditions 
to be satisfied by the flow velocity in the bulk. It turns out that enforcement of the missing condition makes the 
integral equation trivial, in the sense that the only its solution satisfying all boundary conditions is the rectilinear 
front configuration with constant upstream gas velocity. 

Evidently, this fact entails two conclusions. First, the assumptions 1), 2) mentioned above oversimplify the problem, 
and cannot be used to describe nontrivial flame configurations. Second, the limiting procedure proposed by [8| is also 
inadequate. The main purpose of the present paper is to derive the correct equation describing flames in high- velocity 
streams, and to find its nontrivial solutions. For this purpose, the on-shell description of steady flames, developed in 
0, will be used. Its main advantage is that it describes flames in a closed form, i.e., in a form involving only 
quantities defined on the flame front. There is no need to solve the flow equations in the bulk explicitly. In particular, 
it allows one to avoid artificial assumptions about the bulk flow, such as the scaling laws for the gas-velocity and 
pressure in the high- velocity limit, which are adopted in one way or another by the conventional approach. This 
description was initially given for freely propagating flames, but it admits simple and natural extension to anchored 
flames. This generalization is obtained in 

The paper is organized as follows. The Zel'dovich-Scurlock-Tsien approach is critically reviewed in Sec. HT1 which 
starts with a brief account of Cherny's formulation of the high-velocity limit. It is shown, in particular, that this 
formulation reproduces the assumptions 1),2) of Zel'dovich-Scurlock-Tsien approach, and that the only solution of 
the main equation, satisfying boundary conditions, is the trivial solution. The reasons underlying this result are 
identified in Sees. Ill Bt III CI We then go to the on-shell description in Sec. IIIII summarizing the main equations 
derived in and their extension to anchored flames. The reader is referred to for more details concerning 

inclusion of the anchoring system and its analytical description. Sections IIII Al - fill CI are counterparts of Sees. Ill Al 
- Ill CI They discuss the role of pressure, and indicate the place the boundary conditions and evolution equation 
take in our approach. Solution of the on-shell equations is obtained in Sec. IIVI In Sec. IIV Al a large-slope expansion 



of these equations is derived, which extends to curved flames the corresponding expansion constructed in ll| for 
V-flames. Using this expansion, it is shown in Sec. IIV Bl that the main integro-differential equation for the complex 
velocity reduces to ordinary differential equations. Together with the evolution equation, these equations can be 
partially integrated and further reduced to a single second-order differential equation. This is done in Sec. IIV Cl where 
numerical solutions of the derived equations are also found. Section fVl summarizes the results of the work. The paper 
has an appendix which describes in detail transition to the case of vanishingly small anchor dimensions within the 
large-slope expansion. 



II. CRITIQUES OF ZEL'DOVICH-SCURLOCK-TSIEN APPROACH 

A. Cherny's formulation of the high-velocity limit 

We start with brief recalling the results of [H using notation appropriate for the present paper. These results 
provide a rigorous base for Zel'dovich-Scurlock-Tsien approach. Consider a steady two-dimensional combustible ideal 
gas stream in a channel with plane-parallel walls, ignited at a fixed point in the middle of the channel (see Fig. [1] 
where only the right half of the channel is shown). This point (the flame anchor) will be chosen as the origin of 
Cartesian system of coordinates r = (a;, y), with y-axis along the wall, and the initially uniform fresh gas at y — — oo. 
The channel half-width will be taken as a unit of length, while the gas velocity v — (w. u) will be measured in units of 
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the normal flame velocity relative to the fresh gas. Then the flow variables obey the following equations in the bulk 

(1) 
(2) 

(3) 

where p, p are respectively the gas pressure and density. Under the assumption that the flow is essentially subsonic, 
the gas can be considered incompressible. Taking the fresh gas density as a density unit, that of the burnt gas will be 
1/9, where 9 > 1 is the gas expansion coefficient. The flame configuration is assumed to be symmetric with respect 
to the y-axis. More precisely, 

f(x)=f(-x), w(x,y) = -w(-x,y) , u(x, y) = u(-x, y) . (4) 

In particular, the transversal velocity component vanishes at the symmetry axis, w(0, y) — 0, so that this formulation 
applies also to the case of a flame anchored at the channel wall, the point of view taken up in 0] . 

Cherny formulated the high-velocity limit as follows. Denote the fresh gas velocity far from the flame front by U, 
and introduce new (designated with a tilde) coordinates and flow variables according to 

x = x, y = y/U, u~u/U, w = w, p=p/U 2 . (5) 

After that, switch to the new independent variables (y, ip), where ip is the stream function defined by 

dip _ dip 

pu=-, pw = -- B§ . (6) 

Then the limit U — > oo is to be taken under the assumption that the quantities w,u,p as well as their derivatives 
with respect to y,ip remain bounded. Equations d2J) , (J3j) (J6j) thus take the form 

dp ~du 1 dp dx I dx w 

757 = 0, u— + __£=o, —- = — , —: = -. (7) 

dip dy p dy dip pu dy u 

The first two equations give 

~ n 

U T) 

p = p(y), -7r + - = m, (8) 

2 p 

where i{ip) is an arbitrary function. There are also jump conditions to be satisfied across the flame front, which in 
the limit U — ► oo read 

u+ = u- , p+=p-, ip + =ip-, (9) 

where the minus (plus) subscript denotes restriction to the flame front of the flow function defined upstream (down- 
stream). Noting that i(ip) = const upstream (since the incoming flow is uniform), we see that equations ([5]) together 
with the first two equations ^ exactly reproduce the assumptions I) and 2) stated in the introduction. Next, inte- 
grating the third of equations ((7]), and combining its solution with the relations JS]| Cherny arrives at the following 
integral equation for the function y(u—) 

. f uy(u)du (tti — 1) 

"-"/ ^PP : ■ — •- (10) 

which is equivalent to the equations derived by Zel'dovich and Tsien. This completes the proof that the assumptions 
I), 2) of Zel'dovich-Scurlock-Tsien approach are equivalent to the limiting procedure formulated by Cherny. 

But this is not the end of the story. It turns out that the consideration of |8j] is not complete, in that it does 
not take into account boundary conditions for the gas flow. More precisely, these conditions are applied only at the 
end-points of the flame front, but not in the bulk. In terms of the stream function, the condition that the transversal 
component, w, of gas velocity vanish at the walls and at the symmetry axis reads 

■0 = for x = 0, (II) 
ip = l for x = 1, (12) 
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as is seen from equations ([6]) written in the form dip = —pwdy + pudx. Consider the flow upstream. Since i(ip) — const 
there, Eqs. ijHJ) tell us that u is independent of ip, so that integration of the third of Eqs. ([7]) yields 

- * 'X(y), 



where the function X(jj) is to be determined from the boundary conditions pip . (|12|l . Substitution gives 

X{y)=0, u(y) = l. 

Therefore, the flow turns out to be uniform in the whole region upstream of the flame front. Then Eqs. ([8]), ([9]) 
show that so is the flow downstream. Thus, u\ = 1 is the only solution of Eq. (|10[) consistent with the boundary 
conditions at the channel walls. In other words, the limiting procedure proposed by Cherny, and hence the equivalent 
assumptions 1), 2) cannot be used to describe nontrivial flame configurations. 

Technically, this means that the scalings (JSJ) are inappropriate. For instance, the scaling of pressure only seems 
natural. It is suggested by the Bernoulli integral 

- + -(u 2 + w 2 ) — const 
P 2 

which after neglecting w in comparison with large u would give p ~ U 2 . However, it is variable part of pressure that 
only matters. Setting u = U + it, and absorbing U 2 /2 in the "const" on the right gives 

P tt* u 2 

— h Uu H = const . 

P 2 

Now, p ~ U, if u is assumed bounded, or something else if it behaves differently. 



B. Boundary conditions 



Although not always stated explicitly, impermeability of the channel walls (and of the symmetry axis) is assumed in 
the Zel'dovich-Scurlock-Tsien approach, but is not used. As we saw in the preceding section, enforcing this condition 
in the Cherny's formulation makes the flow trivial, leaving the piecewise linear front configuration with constant 
gas velocities up- and downstream as the only possibility. In this formulation, impermeability of the walls and the 
symmetry axis is expressed in the form (fTTj) . (fT2|) . but these conditions are used only at the end-points of the flame 
front, and at the origin. But in fact, things must be just the opposite: transversal gas velocity must vanish at the 
channel walls and the symmetry axis everywhere except a small vicinity of the flame anchor. Indeed, suppose for 
simplicity that the flame is anchored by a cylindrical rod with circular cross-section (as is often the case in practice, 
see Fig. [2]). Then the gas flow is significantly disturbed in a vicinity of the rod: The slowdown of gas elements heading 
the rod leads to appearance of a nonzero transversal velocity component. This happens no matter how small the 
rod radius is. If one follows the trajectory of a fresh-gas element moving near the symmetry axis, its transversal 
velocity rapidly changes near the rod from zero to some finite value at the flame front. This means that in the limit 
of vanishing anchor dimensions (i.e., in the Zel'dovich-Scurlock-Tsien setting), the gas flow is singular at the point of 
its location. 

To be more specific, the presence of the rod in a high-velocity stream can be described mathematically by super- 
imposing the incoming flow velocity v° = (0, U) with the velocity field v d of a dipole located at the origin 

vd= U& { _ 2xyx 2_ y 2 h (13) 

where R <C 1 is the rod radius. Indeed, if the rod is located downstream as is shown in Fig. [2j the flame front bends 
round the rod, and for large U, gets close to its surface, so that the normals to the front and the rod surface coincide. 
To the leading order, the normal gas-velocity is negligible in comparison with U. Hence, velocity of the fresh gas near 
the rod must satisfy 

(»,!/) = 0, (14) 

where v is the normal to the rod. The velocity field v° + v d does satisfy this condition by virtue of the identity 
(v° + v d , r)\ r= R = 0, because r is normal to the rod. On the other hand, if the rod is located upstream (i.e., flame 
is stabilized in the wake of the rod), Eq. (TTJJ is the true boundary condition obeyed by v° + v d identically. 
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The transversal velocity component induced by the rod is large for r ~ R, and rapidly decreases with distance. 
At distances r such that R <C r <C 1, the inner solution describing the flow near the rod is to be matched with 
the large-scale outer solution we are interested in. In particular, matching at the flame front assigns w_ a definite 
value, say wq. Since the transversal velocity in the outer flow is of the order of unity, so is wq. Let us denote by Rq, 
R -C Rq -C 1, the characteristic distance where the two solutions are matched near the flame front. We thus have 
(considering the right branch of the front, x > 0) 

W-\r~Ra = w o , (15) 

where wq is a positive or negative number, while 

w(0,y)=0, y<-R. (16) 
Finally, let us determine the scaling of Ro with U. For r ~ Rq, y = f(x), one has y ~ Ro, x ~ R, and hence 

W-\r~R ~ p4 -K-ftO — d3 • 

Therefore, wo ~ 1 implies i?o ~ U 1 ' 3 R. Taking this into account, we find also 

„<J| ^ P 2 _ rr l/3 

a 

Neglecting [7 1 / 3 in comparison with [7, we conclude that to the leading order, the longitudinal velocity satisfies 

U-\ r ~Ro = U . (17) 

Needless to say that the boundary behavior outlined above cannot be described within the ZePdovich-Scurlock-Tsien 
approach in which transversal velocity component is completely excluded from consideration. 

We have considered the simplest case of a circular rod. In the general case, specific form of the local flow is of 
course different, but the hierarchy of length scales R <C Ro *C 1, as well as relations (fT5")) . (jTTf remain the same. These 
relations will be invoked in Sec. IIII Bl 



C. The role of evolution equation 

The local propagation law of a flame is determined by the so-called evolution equation which gives the normal fresh 
gas velocity, v™, as a function of the flame front curvature and space-time derivatives of the gas velocity at the front. 
For zero-thickness flames, it states that the normal flame velocity is simply equal to the velocity of planar flame. 
Namely, if the flame front position is represented by the curve y — f(x), the evolution equation reads v™ = I, or 

w_ - f'w- = N , (18) 

where the prime denotes ^-differentiation, and N = yl + f' 2 ■ Curiously, this characteristic property which plays 
fundamental role in the whole theory of flame propagation is not invoked in the derivation of the Zcl'dovich-Tsicn 
equation. The reason for this is the already mentioned disregard of the transversal velocity component. The point 
is that the transversal and longitudinal velocity components are strongly coupled by the evolution equation. Indeed, 
it_ and /' are both 0(U), so that the two terms on the left-hand side of Eq. (fT8"|) are of the same order of magnitude. 
Hence, exclusion of w from the list of dynamical variables makes this equation obsolete. In other words, any surface 
of discontinuity with arbitrary propagation law would satisfy Eq. (|10p . provided that its normal velocity is small 
compared to U. 

The only place where Eq. (jT5J) finds application in the Zel'dovich-Scurlock-Tsien approach is the relation between 
the flame front length and the incoming stream velocity 

U = f(l), (19) 

which follows from the mass conservation condition 

l l 

U = J dxy/T+f*v1 « J dxf = /(I) - /(0) . 



Of course, this relation is independent of the particular dynamical model employed, and is valid in the high-velocity 
limit whatever the structure of the incoming flow. In the Cherny's formulation, Eq. (|I8[) takes the form dip_/dy = 1 , 
and after integration under conditions (fTTj) . (jX3]) yields Eq. (II9p . These equations are not used in the derivation of 
Eq. (USD- 
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III. ON-SHELL DESCRIPTION OF FLAME PROPAGATION 

As shown in 0, [T3], the system of bulk flow equations and jump conditions at the flame front can be reduced to 
a single complex integro-differential equation relating values of the flow variables at the flame front (their on-shell 
values), so that explicit solving of the bulk equations (which is the most difficult part of consideration) turns out to 
be unnecessary. This equation reads 



2( W _)'+(l + iM){[ W r-^^} = 0. 



(20) 



Here, u = u + iw is the complex velocity, [u>] = lu + — w_ its jump across the front, v™ is the normal velocity of burnt 
gas at the front, a+ is the on-shell value of vorticity produced by the curved flame, and the operator j£ is defined on 
2-periodic functions by 

+i 

fc») (x)= 1 + i f' (x) jdrj a( V ) cot {^( V - x + i[f ( V ) - f(x)})} , (21) 

-l 

the slash denoting the principal value of integral. It satisfies the important identity 

3< 2 = -l. (22) 

The functions er+,i>? as well as the velocity jumps at the front, entering equation (|20p , are all known functionals of 
on-shell fresh gas velocity [see, e.g., [H, For zero-thickness flames, 

6? = », M = ^. M = "/'^. (23) 

= -W<«- + "'-)'■ < 24 » 

Together with the evolution equation (|18p. Eq. (j2"D|) constitutes a closed system of three equations for the three 
unknown functions W-{x),U—{x), and fix). 

Equation ([20]) describes freely propagating flames, but can be easily modified to take into account the presence of 
the rod. This generalization is given in [Tlj . and reads 

2 („_)' + (l + ift) { M ' - ^>±^±} = 2 (.1)' , (25) 

where co d is the complex velocity of the dipole ([T3]) . Since this field satisfies the symmetry relations (j4|), the boundary 
condition (|16[) is still met. We assume for definiteness that the rod is located in the downstream region (as is usually 
the case in practice; this assumption is inconsequential, Cf. footnote [T]) . Then uj- and u>l satisfy 



(l-itt) (w_)' = 0. (26) 
(l (wl)' = 0, (27) 



which express analyticity and boundedness of the functions u>(z), ui d (z), z = x + iy, in the upstream region. Equations 
(|26|) . ((27)) are consistent with Eq. ([25]) by virtue of the identity ([22]) . 

We will now consider the problem of anchored flame propagation within the on-shell description in the light of the 
issues discussed in Sec. Ql] 



A. The role of gas pressure 



It was already mentioned in the Introduction that specifics of the two-dimensional case make the issue of flame flow 
nonlocality especially nontrivial. The region near the rod is characterized by large velocity gradients, and therefore 
affects significantly the pressure field far apart from the rod. Indeed, pressure is determined by the Poisson equation 

Ap = -p(V(vV)v) , 
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while the Green function of the Laplacian is proportional to In r. The flow near the anchor thus can be expected to 
have a strong nonlocal influence on the flame structure. It is one of the advantages of the on-shell formulation that 
it reveals the completely subordinate role pressure plays in determining the front structure. In fact, the variable p 
does not appear in Eq. (|20[) . In particular, no assumption such as that contained in the point 1) of the ZePdovich- 
Scurlock-Tsien approach, or the scaling prescription ([5]) of the Cherny's limiting procedure, is needed in the on-shell 
description. Also, exclusion of pressure from consideration weakens the above argument concerning nonlocality of the 
anchor impact. 

B. Strategy of solving Eq. (|25|) . Boundary conditions 

We are concerned with the situation where the rod radius, R, is much smaller than the channel width, and interested 
in the large-scale front structure, i.e., its structure at distances large compared to R (the outer solution). We observe, 
first of all, that the presence of the rod is described simply by adding the dipole field to Eq. (|20p. which is noticeable 
only in a small vicinity of the origin. This fact indicates that despite the pressure argument given in the preceding 
section, the detailed flow structure near the rod may be actually unimportant. This would mean that the rod can be 
considered point-like, thus greatly simplifying account of its influence on the global flame structure. Considerations of 
the subsequent sections will show that this is indeed so in the high- velocity limit. Namely, it turns out that Eq. (|2"0|) 
reduces in this limit to an ordinary differential equation. This naturally opens the following way of solving Eq. (I25|) : 
we consider this equation at x's such that r > Rq, and extract the leading terms of the high-velocity expansion. 
Detailed analysis reveals that for such x's, the dipole contribution is negligible, so that Eq. P5|) reduces to Eq. (|20p . 
Exclusion of the region near the rod implies that this equation must be supplemented by an auxiliary condition at 
a point xo corresponding to the matching region r ~ R a (see Fig. [2]) . Since Eq. (|20p is complex, and involves only 
quantities defined at the flame front, it requires two real auxiliary conditions expressed in terms of on-shell quantities. 
As such, we take the relations (fT5)) , (fTT)) in which the numbers Wq,U play the role of parameters specifying given 
large-scale solution. In the limit R — ► 0, these relations take the form 

w-{0 + )=w , (28) 
u_(0+) = J7, (29) 

with the understanding that the functions W-(x) = w(x, f(x)), U-{x) = u(x, f(x)) are solutions of Eq. (f2"0"|) considered 
on the semi-open interval x £ (0,1], and the quantities w~(0 + ), u_(0 + ) are their right limiting values for x — » 0. 
Furthermore, taking the limit R — > leads obviously to the following condition for the function f(x) 

/(0) = 0. (30) 

Finally, there is also the usual condition of vanishing of w at the wall x = +1, 

iu_(+l)=0. (31) 

If a solution of Eqs. (JTSJ, (f2"0")) satisfying the above conditions is found for x £ (0, 1], then the rule (HJ) gives it on the 
semi-open interval x £ [—1,0) as 

f(x) = f(-x), w-(x) = -w-(-x) , u-(x) = u-(-x) , (32) 

in particular, w—(x) satisfies u>_(0~) = — wq, so that this function has a discontinuity at x = 0. 

In connection with the above procedure of finding the large-scale solution, the following circumstance should be 
emphasized. We are able to formulate this procedure in a closed form, without the need to construct the local solution 
near the rod explicitly, just because Eq. (f2"0")) relates only functions defined at the flame front. This would not be 
possible in the conventional approach based on explicit solving the bulk equations. Indeed, as we saw in Sec. Ill Bl 
the transversal velocity of gas elements moving near the rod rapidly changes from zero to some finite value at the 
flame front. Thus, if we were to construct a bulk outer solution satisfying the matching and boundary conditions (|15p . 
(1161) . this solution ought to describe the rapid change of w from zero to wq ~ 1 over a distance ~ Rq. One might try 
to avoid this complication by shifting the matching region farther from the rod, but then the dipole field would be 
completely neglected. In effect, the anchor disappears from the description, and the solution becomes trivial. This is 
exactly what happens in the Zel'dovich-Scurlock-Tsien approach. In our approach, on the contrary, only longitudinal 
velocity induced by the dipole is negligible in the matching region, but not the transversal. 

At last, it is worth noting that the parameters wq,U specifying the large-scale solution are independent of each 
other. This fact becomes evident when we note that instead of the rod with circular cross-section we might use a more 
complicated cylindrical shape. Then for the same value of U, matching of the inner and outer solutions would give 
a different value for wq. We conclude that the large-scale solutions form a two-parameter family. This is in contrast 
with the Zel'dovich-Scurlock-Tsien approach where the single parameter U completely determines the solution. 
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C. Evolution equation 

Having formulated the outer problem in a closed form, we may now note that the experimentally observed anchored 
flames usually have fairly smooth, highly elongated front configurations. Hence, neglecting the small regions near the 
rod and the channel walls where finite-thickness effects are only important, we can use the evolution equation in the 
simplest form (TT8|) applicable to zero-thickness flames. Restricting ourselves again to the right half of the channel 
where the front slope is positive, and taking into account that /' = 0(U) [Cf. Eq. (fTO|) ] this equation can be rewritten 
as 

u-=f'(w- + l), x>0, (33) 

the correction term being of the relative order 0(1/U 2 ). 

It should be emphasized that neglecting the regions near the rod and the channel walls makes boundary conditions 
for the function f'(x) itself unnecessary. Such conditions are only needed in establishing the detailed structure of 
these regions, which is directly related to the fact that the finite front-thickness effects determining this structure 
are described by equations of higher differential order. This means that the slope of the function f(x) describing the 
outer solution is to be considered large everywhere on the semi-open intervals x G (0, 1] and x G [—1,0). The points 
x = ±1 are included here by continuity, but not the point x — which represents the true singularity of the flow, 
where /' is undefined. 

IV. SOLUTION OF EQ. ([20]) IN THE HIGH- VELOCITY LIMIT 

A. Large-slope expansion of the IK-operator 

The nature of nonlocality of steady flames is encoded entirely in the structure of the operator JC, and of major 
importance is the fact that this operator greatly simplifies in the high- velocity limit. For large U, the front slope is 
also large, |/'| ~ U, so the argument of cotangent in Eq. ([?!]) has a large imaginary part for almost all values of the 
integration variable. Therefore, one can write 

cot {| (77 - x + i[f{r,) - f{x)])} w -»x(|»7| - \x\) , (34) 

where xi x ) is the sign function, 




+1, x > 0, 
-1, x < 0. 



This approximation is valid for all r\ except two small regions near r\ — ±|x| . More precisely, taking into account that, 
for real ai^, 

e (o 2 -iai) g— (a 2 -iai) , . 

COt(ai + ia 2 ) = —i —, : — ; ; : — c = -ix( a 2) + O ( e _2|a2 ) , 

we see that Eq. (|34[) holds true, with an exponential accuracy, everywhere except 

ij:|» 7 |e(|a;|-a,|a;| + <J) ) 

where 5 = 0{l/U). 

To develop an asymptotic large-slope expansion of 3i, let us choose a real e > such that 

e < 1 , Ue » 1 . (35) 
Then the integral in Eq. pTj) can be rewritten, for x > 0, as 

+i 

drj a{rj) cot (rj - x + i[f{rj) - f(x)])} 



—i 



drj a(r)) cot || (77 - x + i[f(rf) - f(x)})} . (36) 



+ 1" 



-x+e 



X + E 



drj a(rj)x(\v\ - x ) 
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Notice that in the last term on the right hand side of Eq. (|36|) , only one of the two integrals is defined in the principal 
value sense. As such, it is proportional to a'(x). It is not difficult to see that contributions of this kind give rise to 
terms of the order l/U 2 . This is because expanding the function a(rj) around x brings in an extra small factor (rj — x). 
Below, we will need "K expanded up to 0(l)-terms, so the principal-sense integral can be neglected. The other integral 
can be evaluated as follows, using continuity of a(rj) 

— x+e +e 

(IT ~\ f ( 7T f f (^c) 

dr] a(r]) cot |- (t] - x + i\j(r\) - f(x)])j = -ia{-x) / dy coth j — - — y + Trix 

—e 

o +7rf'(x)e/2-\-7vix 



-ia(-x) In shy 

nf'{x) 



-7T /' (x)e/2-\-7rix 



In view of (|35[) , In shy can be replaced, with the exponential accuracy, by y and — y at the upper and lower integration 
limits, respectively. Taking into account also that for x — » + , arg(shy) gains — 7r, we find 



+7r/ y (x)s/2-\-7rix 

= Tri{2x - 1) . 

— irf r (x)e/2-\-Trix 



In shy 

On the other hand, replacing the cotangent in the last term in Eq. ()36[) by the sign function gives zero within the 
same accuracy 

-x+e +e 

d V a(v)x(\v\ -x) = a{-x) j dij x(v) = • 

—x—e —e 

Using these formulas in Eq. (|36p . and then putting it in Eq. (|2~Tj) gives finally 



($Ca) ( 



,., = ,/•',,.). / d V a ^ )+ < ^ x{v -\ x \) +ia (-x)(2\x\-l)+o(^ 



(37) 



This result is written in the form applicable to negative as well as positive x, which can be verified by noting that 
V}i is invariant under the combined operation of coordinate inversion (x — > —x) and complex conjugation, as is seen 
from Eq. pip . In particular, the asymptotic action of IH on the derivative of a function a{x) continuous at the origin 
and satisfying a(— 1) = a(+l), is 

(%a') (x) = (f'(x) - i) {a(-\x\) - a(\x\)} + ia'(-x)(2\x\ - 1) + O ^ , (38) 

where the prime now denotes derivative of the function with respect to its argument, a'{y) — da(y)/dy. It turns out, 
however, that the formula (|38|) remains valid even if the conditions a(CP) = a(0 + ), a(— 1) = a(+l) are not met. In 
particular, a(x) can be discontinuous at x = 0, so that its derivative is singular at the origin. This important fact is 
proved in the appendix. 

Formula (l38|) was derived for x's where the front slope is large, i.e., for all x except small regions near the rod and 
the channel walls, where front curvature is large. Neglecting these regions as we did before, we can say that Eq. (|3"5|) 
is valid on the semi-open intervals x G (0, +1] and x £ [—1,0) (see Sec. IIII C|) . 

The following comments concerning the structure of Eqs. (|37| . (|38j) will be useful in subsequent applications. First, 
it is seen that the result of the action of !K depends essentially on parity properties of the function a(x), namely, 
3ia = 0(U), if a(x) is even, and 3ia = 0(1), if it is odd. Second, it should be noted that although the identity 
3i 2 — — 1 is valid whatever the shape of the flame-front, in particular, in the large-C/ limit, it cannot be verified using 
the expression on the right of Eq. (|37p . already because the composition of its leading term with the undetermined 
remainder 0(U) o 0(1/U) = 0(1). 
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B. Reduction to the system of ordinary differential equations 

1. Equation for the transversal velocity component 

We now go over to proving the results announced in Sec. IIIIBI First off all, let us determine the orders of various 
terms in Eq. (|2"5|) within the high-velocity expansion. As we know, u_ = 0(U), f = 0(U), so Eq. (|3"3")l tells us that 
W- = 0(1). Using these estimates in the expressions (|23|) . (|24|) shows that 

[it] =0(1/U), [w]=0(l), u + = 0(U), w+ = 0(l), v+ = 0(U), Na+ = 0(U 2 ). 

Therefore, one has for the braces in Eq. (|25[l 



Re { [u]' - ^±^±) = (U) , Im ( K - M^±\ = 0(1) . 



Since Re{...} is an odd function of x, while Im{...} is even, formula (|37[) shows that the leading term in the real part 
of %%{...} is 0(U) . Thus, the real part of the left hand side of Eq. (f^5|) is 0(U). At the same time, as we saw in 
Sec. (jll B[) . the real part of the dipole velocity field is 0(1) in the matching region, and negligible far apart from the 
origin. Therefore, to the leading order of the high-velocity expansion, the real contribution to the right hand side of 
Eq. (|25p can be omitted. However, things are quite different for the imaginary contribution. In this case, the above 
estimates and formula l|3"Tj) show that both sides of Eq. ([25]) are 0(1). Moreover, expansion of !H only up to 0(1)- 
terms is actually insufficient for the purpose of extracting the imaginary part. Indeed, the undetermined remainder 
in Eq. (|37[) is 0(1/U), and may contain real as well as imaginary parts. Since the argument of the IK-operator is 
0(U), this remainder gives rise to terms of the order 0(1/U) -0(U) = 0(1). These two complications can be overcome 
by resorting to Eq. (|26p which expresses potentiality of the upstream flow, and can be considered as a consistency 
condition for Eq. ([25]) . Repeating literally the above arguments, 1 one sees that the real part of Eq. (|26|) can be 
consistently extracted with the help of the expansion obtained in the previous section. Namely, using the formula 
(|57)> yields 

u'_(x) - f(x){w-(\x\) -w-(-\x\)} + u'_(-x)(2\x\ - 1) = 0, 
or since u'_(—x) = —u'_(x), W-(—x) — —w(x), 

u' (x) 

w _(\x\) = (l-\x\)-±-L. (39) 
f(x) 

By the construction, this equation is valid for x S [—1,0) U (0, 1]. In particular, we see that the boundary condition 
(f3"Tj) is satisfied automatically. 

Thus, we proved that in order to find the large-scale solutions of Eq. ([23)1 , it is sufficient to consider Eq. (|20|) . 

2. Equation for the longitudinal velocity component 

Turning back to extracting the real part of Eq. (|20p , we have to consider the question concerning the contribution 
of the small region near the rod. As we saw in the preceding section, the dipole field on the right of Eq. (|25p can 
be neglected for i^O. However, this does not settle the question, because integration on the left hand side extends 
over all x including zero. The first term in the braces is a derivative, so that we can use Eq. (|38p to find how it 
is transformed by "H. Note that since this term is the derivative of a function discontinuous at x = 0, it contains 
contribution proportional to the Dirac (^-function. Indeed, one has from Eq. (|23p 



1 .ifwt-Ln-w-iMx), 



N N 



1 For a flame stabilized in the wake, Eq. 1261 1 is replaced by (1 — i6-Cj )' = 2 (u'L)' 11] • This makes no difference for the present 
analysis, as the dipole field is negligible on the same grounds as in Eq. 1251 1. 
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so that 

[w]' = -2i{9-l)S(x) . 

However, it is proved in the Appendix that the formula (|38p is still applicable in this case. As to the second term in 
the braces, it also contains a derivative factor, viz., (u_ + w_)' coming from cr + , but this time this is a derivative of 
an even function. However, it would be premature to conclude that this term does not contain a ^-contribution. The 
point is that (ut + w 2 .)' is multiplied by w+ whose imaginary part is an odd function. Let us trace the development 
of the quantity Q = Nw + cr + /v'] r near the rod (see Fig. [3]). At the matching point x = —xq on the left of the rod, trL 
is large but changes slowly, while w + = —wo + (6 — 1) = 0(1). Also, if the bulk transversal velocity of fresh gas is not 
too large, and \v+\ increases away from the anchor (pressure normally drops down along the stream), then w + > 0, 
{v\)' < 0, and so 

< Q = -™+^7^ (!*«+)' = 0(1) • (40) 

<T+ 1 increases near the rod as the result of gas slowdown caused by the rod, and for |x| < R, Q becomes 0(U/R). 
Here, the front curvature is large, and the zero-front-thickness expression (|40j) can be used only for a rough estimate. 
It shows that Q is negative in this region, because u>+ < 0, and er+ > 0. Q rapidly turns into zero at x = 0, because 
both w+ and ct + vanish at the origin. For positive x < R, Q is again a negative 0(U/R) quantity, since w + > 0, 
er + < 0. At larger x's its modulus decreases, and Q becomes 0(1) at the matching point x — xq on the right of the 
rod. From the large-scale point of view, this behavior means that Q contains a term qd(x), with a negative coefficient 
q. The exact value of q can be found, of course, only if the inner solution is known. We arrive at the conclusion that 
the expression in the braces in Eq. (|20p can be written as 

Au 2 + w 2 Yuj+ 

-iqS(x) + (9-l) K ~ 2v ~> + , 

where q = q + 2(0 — 1), and it is understood that the (5-contribution is excluded from the second term. In other words, 
this term is calculated using the functions u_,w + etc. that describe the outer solution. Taking into account that 
v\ = v 2 _ + 8 2 — 1, using Eqs. (|57|) . (|38[) . and extracting the real part of Eq. (f!20")) gives, to the leading order, 

l 

u>_(x)(l + a\x\)-p>(x)-^f(x) jdr,^-[w^r,)-a}x(v-\x\) = 0, a = 0-1. 

o 

This equation involves the unknown parameter q. To get rid of it, we divide the equation by /', and then differentiate 
it with respect to x. The result is the following ordinary differential equation 



d 
dx 



(1 + a\x\) 



fix) 



u' (x) 

■ a^^[w-(x) - a] = . (41) 

u_ (X) 



Together with Eqs. (|33|) , ((39)) it constitutes the system of three ordinary differential equations for the three functions 
U-(x), W- (x), f(x). Evidently, it requires three initial conditions which are Eqs. ^ZE\i - (|30|) . 



C. Reduction to a single differential equation. Numerical solutions 

Introducing an auxiliary function 

u'_ dii- 

the system ([33)) . (|3l?)) . (|H)) can be rewritten as an ordinary differential equation for <p(x): 

d r , (1 — x)<x> — a , 

Tx^ 1 + a ^ + ^ \l-x)l + l =^ X> °- (42) 

The initial condition for <p follows from Eqs. ([28)) . ()39| : <f(0 + ) = wq. We note also that the functions U-(x), f(x) are 
related by a simple algebraic equation. One has from Eqs. (|33)) . ([39)) 



(1 - x)u'_ + /' , (43) 
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[(1 -!)«_]' + /' = 0. 
Integrating this equation, and using the initial conditions ([29)) , (|30|) gives 

/ = Z7 - (1 - x)u- . (44) 
Finally, combining Eqs. (|43|). (|44l one can express <£> in terms of / 

£7 - /(x) 1 

^ ' (1 - X) 2 /'(X) 1 - X ' 

Substitution of this expression into Eq. (|4"2")) leads to a second-order differential equation for the front position. The 
corresponding initial conditions follow from Eqs. lj2"8)) - ([50)1 . and 

/(Q)=Q, /'(Q+) = -^ T . (45) 

Numerical solutions of these equations for various values of U, wq and 6 are plotted in Figs. 0] - [HJ They have 
the following general features. First of all, the function f(x) is monotonic in all cases (in each of the two channel 
halves), as was assumed throughout our consideration. Second, flames in which the fresh-gas flow diverges near the 
rod (wq > 0) are convex towards the incoming flow, while those with convergent fresh-gas flow (wq < 0) are concave 
(recall that we deal here with the large-scale solutions, characterized by distances r > i?o from the rod; for |x| < R, of 
course, wq is positive in any case). Numerical analysis shows that in the latter case, solutions with the given negative 
wo exist only for sufficiently small values of the gas expansion coefficient. For example, solutions with the fairly small 
value wo = —0.1, one of which is shown in Fig. [6) disappear at 6 « 3.5, and this threshold value is independent of U. 
Finally, solutions with wq > are characterized by monotonic increase of w_ as one moves from the rod to the wall. 
The overall velocity rise is substantial - typically two to four times. This is what normally observed in experiments. 
Solutions with wo < are anomalous in this respect, as U- slightly decreases away from the rod. They are most likely 
unstable. 

At last, solutions with wo = are trivial. Indeed, in this case ip(0 + ) = 0, and Eq. (|42f tells us that also f'(0 + ) = 0. 
Repeated differentiation of this equation then shows that all higher derivatives of tp also vanish, i.e., ip{x) = 0. Hence, 
it_ = const = U, and Eq. [(li)) gives f(x) = Ux, x > 0. In the light of the discussion given in Sec. Ill Bt it is natural 
that the case wq = reproduces the result of Sec. Ill Al 



V. CONCLUSIONS 



The results obtained in this paper provide consistent description of steady anchored flames in high-velocity gas 
streams in channels. Given the values of the incoming flow velocity and its transversal component near the anchor, 
the formulas derived in Sec. IIV CI allow simple determination of the flame front shape and on-shell gas-velocity. A 
practically more convenient may be "geometrical" parametrization using the ordinates of the front end-points and its 
slope at the origin, which is related to the initial one by Eqs. (TIT))) , (|45[) . 

The remarkable fact revealed by the above investigation is that the flame structure in a high- velocity gas flow obeys 
an ordinary differential equation. In other words, this structure turns out to be local in the usual sense: behavior 
of the flame front slope and gas velocity in an infinitesimal vicinity of a given point is determined by their values 
at this point. The sole role of the anchor is to provide an initial condition. This result answers the question as to 
the nature of nonlocality of the anchor influence on the flame structure: Although detailed structure of the flame 
holder is immaterial for the properties of the large-scale flow, the flow distortion it causes ultimately determines the 
whole flame configuration. We saw in Sec. IIIBI that mathematically, the presence of the anchor with vanishingly 
small dimensions signifies existence of a singularity in the bulk flow solution. The failure to recognize this fact is 
what makes it impossible to consistently describe nontrivial flame configurations within the ZePdovich-Scurlock-Tsien 
approach. Actually, this defect is inherent to this approach as it neglects the transversal gas velocity, while the role 
of this component is crucial in describing the anchor impact. 

Finally, the role of vorticity in the formation of curved flames is to be emphasized. It is described by the second 
term in the braces in Eq. (|20[) . while the first term (the complex velocity jump) corresponds to a purely potential 
contribution. As we have seen in Sec. IIVB 21 the latter falls off from the equations describing the large-scale flame 
structure. One can say that formation of the steady flame pattern in a high-velocity stream is governed by the 
vorticity generated in the curved flame front. Therefore, it cannot be described within potential-flow models such as 
suggested in [l4j]- 
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APPENDIX A: EXTENSION OF EQ. ([38]) TO DISCONTINUOUS FUNCTIONS 



If the function a(x) in Eq. (|38|) does not satisfy conditions 

a(0+) = a(0 _ ), a(+l)=o(-l). 



(Al) 



its derivative is singular at x = 0, ±1, and the integration by parts used in the transition from Eq. (|37j) to Eq. (|38j) 
is ambiguous. We recall that the functions describing the true flame configuration are actually smooth and periodic, 
and hence satisfy the conditions (|A1[) , whereas discontinuities arise as the result of simplified description. Therefore, 
in order to correctly evaluate the integral, one has to turn back to the exact formula (f2"Tj) in which all the functions 
involved are smooth, and apply it to a function A(x) satisfying (|A1|) . whose behavior near the rod or channel walls 
looks discontinuous from the large-scale point of view. More precisely, A(x) is supposed to vary rapidly for |x| < R <C 1 
and near the walls, but normally at the intervals xo < x < 1 — xq and — 1 + xo < x < —x$, where it coincides with 
a(x). Here the positive numbers R, x are such that R < Xq <C 1; they have the same meaning as in Sec. IIIIBI Thus, 

lim A(x) = a(x) . 

(R — ► implies that xo also goes to zero.) Neglecting the anchor dimensions means that the action of "K on a' is 
defined as 



(tia') (x)= Um {(jL4')}( a 



(A2) 



To find out how % acts on the derivative of A(x), we replace a by A in Eq. ([37|) . and integrate the right hand side by 
parts 



"HA') (x) = 



l + if'(x) 



dr, A'( V ) cot {|(r? ~x + i[f(r,) - f(x)])} 



1 d 

2dx 



dr, [1 + if(v)}A(ri) cot {^(r, -x + i[f{rf) - /(x)])} 



(A3) 



The boundary terms vanish here because the integral kernel is 2-periodic, and A(x) satisfies A(— 1) = A(+l), by the 
assumption. The function a{x) is allowed to have only a finite jump, and so is the slope, /', of the limiting form of 
the front. Therefore, the last integral in Eq. (|A3|) . in which all functions are replaced by their limiting expressions, is 
well-defined, representing a continuously differentiable function for all |x| G (0, 1). Thus, we can write 



1 d 



+i 



fim { (?U') } (x) = - — j dr, [1 + if'ivMn) cot {^(r, - x + *[/(„) - /(*)])} , 



it being understood that / in the integrand is used in its limiting form. 

Next, we go over to the large-slope limit. The right hand side of the last equation can be evaluated in this case in 
exactly the same way as we arrived to Eq. (|37j) . Comparison with Eq. (|36|) shows that the role of the function a(j]) 
in this equation is now played by [1 + if (rj)]a(ri), the only difference being that the large factor /' comes from the 
integrand, rather than from the pre-integral factor in Eq. (|2 1 [) . Taking this into account, we readily find 



(Ma) (a 



1 d 

2dx 



— ia(— x)(2|x| — 1) 



dr) {a(r))[f W + a{-r))[f' (~v) - i]}x(v ~ M) 



= -/'(|x|) X (x) {a(|x|) - a(-|x|)} + i X (x) {a(|x|) + a(-|x|)} 



-2ia(-x)x(x) + ia'(-x)(2|x| - 1) 
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Using the obvious identities f'(\x\)x(x) — f'(x), x(x){a(|x|) +a(— |x|) — 2a(— x)} — a(\x\) — a(— |x|), we finally obtain 

(&a') (x) = (f'(x) - i) {a(-\x\) - a(\x\)} + ia'(-x)(2\x\ - 1) , 

which is exactly Eq. (|38[) . as was to be proved. We also observe that the result is independent of the particular choice 
ofA(x). 

Moreover, it turns out that this formula is valid not only on the open intervals x G (—1,0) U (0, 1), but in the 
whole channel domain x G [—1, +1], if the derivatives of functions discontinuous at x = are understood in the sense 
of distributions. Having in mind possible future applications, let us prove this fact. Note first of all, that if a(x) is 
discontinuous at x = 0, i.e., a(0 + ) — a(0~) = [a]o ^ 0, then for the function b(x) — a(x) — [a]ox( 2; )/2, one has [b] = 0, 
so that Eq. ([38]) is valid for b(x). Writing a(x) = b(x) + [a]ox(x)/2, we see that since !H is a linear operator, it is 
sufficient to prove the above statement only for the sign function. Let X(x) be its smooth approximation. Take a 
test function <fi(x), i.e., a smooth function that slowly varies for x ~ R, and integrate it with Eq. (IA3j) over interval 
-A < x < +A, where A is such that R < A < 1. We get 

A +1 



/ 



dx<t>{x) (MX 1 ) (x) = 1 jdr, [1 + if{ri)]X{r,) 
-l 

x 0(A) cot {|(r? - A + $[/(»?) - /(A)])} - 0(-A) cot {|(, ? + A + i[f( V ) - /(A)])}] 

A +1 

-\ J dx<t>\x)j : dr ] [l+if(T 1 )}X(r 1 )cot{^(r 1 -x + ilf(r,)-f(x)})}. 

-A -1 

All integrals on the right are well-defined in the limit R — > 0, so that X(rj) can be replaced by x(rj). Then the 
^-integrations are readily done because the primitives are lnsin{-}. For example, 

+i 

j dr, [1 + z/'W]^) cot {|(ry - A + »[/(,,) - /(A)])} 

^ - J dlnsin - A + i[/(r?) - /(A)])} = 2[/(l) - i] - 4[/(A) - iA] 

.0 -1 . 

where it is taken into account that the front slope is large for \q\ ^> R. A simple calculation gives 

A 

dx<t>{x) [viX'^j (x) = -2i0(O) - 2 / dx0(x)[/'(a;) - i] . 



-A -A 



Finally, since 4>(x) is independent of R, the limit of this equation for R — > can be written using the definitions of 
Dirac (5-function and (I A2D as 



y dx0(x) ^x') (a:) = -2i J dxS{x)<j){x) - 2 y dx</»(x)[/'(x) - i] , 

-A -A -A 

which in view of arbitrariness of cf>(x) yields 

(% X ') {x) = -2i5{x)-2[f {x)-i\. 

By virtue of the relations x'( x ) = 25(x), |x|<5(x) = 0, understood in the sense of distributions, this is just Eq. (|38|) for 
a = X- 
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